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Abstract 

Peierls gap is analyzed in case of a two-dimensional lattice under the 
influence of a magnetic field, in a tight-binding approximation. By using 
a non-analytic class of potentials, such as the Kohmoto potential in the 
Harper model, splitting effect occurs in the energy bands. In the metal- 
insulator transition point, the charge carriers become separated due to 
their energy, releasing and expanding the Peierls gap. As a result, the 
energy bands around the Fermi level become localized in case of the 
electrons and delocalized coresponding to the holes, since their energy 
become lowered. These facts are supported by numerical investigations. 
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1 INTRODUCTION 

Studies regarding two-dimensional electron systems under the simultane- 
ous influences of a periodic potential and of an external homogeneous magnetic 
field have received much interest in the last four decades [H El [H E] • This 
kind of systems are important in solid state studies, because the nowadays 
nanoscale technologies are based on the band engineering concept O [7] and 
on the thermodynamic formalism [8j. In order to build heterostructures with 
a precise distribution of the energy bands, we have to understand the way in 
which they appear. 

The energy bands structure of the carriers in magnetic fields was pre- 
dicted in 1964 by Azbel, while in 1969 Langbein [S] analyzing the Landau 
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levels, discovered "the butterfly pattern" , in fact a new and profound connec- 
tion between the admissible energy states and the form and size of the latticed 
pierced by the magnetic field. This phenomenon was analyzed numerically by 
Hofstadter [TO] in 1976. He obtained the self-similar structure of the energy 
bands, with a shape of a butterfly. 

It has been found that electrons moving under the simultaneous influence 
of a periodic potential and of a magnetic field display a series of mini-bands, 
separated by corresponding mini-gaps. The energy spectrum due to its unique 
shape, represent even today a remarkable achievement of the theoretical in- 
vestigation in physics. 

The interpretation given by Hofstadter to the numerical results, high- 
lighted the fact that for a given value of the anisothropy parameter, the system 
is at a crossroad between the conductive and insulator states. This self-dual 
point between localized and extended states, have been described by Andre 
and Aubry [11] too. In this critical point appears a gap described by Peierls in 
1955 [12] known today as the Peierls gap |13} I14j. This represents ever since 
a cornerstone to any studies connected to transitions between conductive and 
insulator states [15] . 

The aim of this paper is to investigate numerically the energy bands dis- 
tribution of a system in which the Peierls gap create a separation between the 
charge carriers, leading from conductive to insulator states. For this purpose 
we intend to analyze the Harper model with the Kohmoto quasi-periodic po- 
tential. Similar studies have been done by Hiramoto and Kohmoto |16j by 
resorting to a generalized Harper model. 

The paper is organized as follows: After the Introduction, in Section II 
we present as a short review of potentials of a two-dimensional lattice. In 
the following one shows the Harper equation with the Kohmoto quasi-periodic 
potential. The numerical results are discussed in Section III. In the last Section 
we present the Conclusions of our research. 

2 Model description and numerical results 

In order to analyze the Harper equation with the Kohmoto quasi-periodic 
potential, we have to reiterate the concept of potentials. That's important 
because we are interested to identify the way in which carriers cross the envi- 
ronment from a point to the other, depending on their conductive properties. 

Let us start from a one-dimensional system in a tight-binding approxi- 
mation, submitted to an external and homogeneous magnetic field [17\ I18|. It 
can be viewed as a conductive media of a free electron gas. By generalizing 
the above chain into a two-dimensional square lattice, we can describe any 
particular direction using the Miller indices [u v w] [19j . In this way, three 
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kind of periodic potentials arise. 

The first kind concerns the periodic potentials [20] which corresponds 
to the directions [0 1 0] parallel to the x-axis, [1 0] parallel to the y-axis 
and [f f 0] for the first and second bisectrices. These directions work in 
conjurrection with the NN (nearest neighbor) and the NNN (next nearest 
neighbor) hoppings. The periodic potentials are described in terms of the 
Bloch functions. They have a rational number as slope, but not localized 
states. The second kind concerns the super-periodic potentials [21], which 
occur to large u of the corresponding Miller indices [u 1 0] . These potentials 
have a large integer as slope and lead to localized states. 

The third one, the quasi-periodic potentials [22J [23] have a regular but 
never exactly repeating pattern with an irrational number as slope. In this 
later case, the translation invariance is lost. One deal instead with local rota- 
tional symmetries. It is a kind of long range order potentials, with similarities 
between tiny bits appearing in different places but slightly modified, from one 
to the other. The quasi-periodic potentials can be found for instance, between 
[u\ 1 0] and [u2 1 0] directions. 

At the beginning these kind of potentials where only strange entities 
without any applicability. Howewer, between 1962 and 1984 they became the 
major description tool of a sort of quasi-crystals [21] ■ One of the most in- 
triguing properties of these solid state compounds are connected with their 
behavior at low temperatures. It is known that any conductive medium re- 
sponds positively to the influence of low temperatures [25] by increasing its 
conductivity and giving birth to the superconductor state. The quasi-periodic 
crystals behave vice-versa, by transforming from a conductor state into an 
insulator one. Such a crystal is provided by synthetic organic materials like 
TTF - TCNQ (Tetra Thia Fulvalene - Tetra CyaNoQuinodimethane) [26], 
which have been discovered between 1962 and 1973. 

Besides the organic compounds, there are other kind of crystals which ex- 
hibit a forbidden five-fold symmetry, such as the Al-Pd-Mn icosahedral quasi- 
crystals [27], discovered in 1984. This alloy, formed with 3 metals, each of 
them with good conductivity, is in fact almost an insulator. 

An experimental attempt to understand the way in which electrons flow 
through a quasi-periodic crystal have been done by Torres et al. [28] in 
2003, resorting to a pulse propagation with a shallow fluid covering a quasi- 
periodically drilled bottom. In order to implement the quasi-periodicity, we 
have to resort to the Harper equation [29J, given by 

VWl + V'n-l + VnlPn+1 = E ^n i l ) 

where ip n is the wave functions at site n. Here V n is a periodic function with 
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the role of potential 

V n = 2A cos 2iT(n(3 + , (2) 

A being the anisothropy parameter and E the allowed energies. 

In other words we deal with a two-dimensional system under the influence 
of an external and perpendicular magnetic field. The quotient (5 = P/Q = 
4>/<fto is the commensurability parameter expressing the number of flux quanta 
per unit cell, P and Q are two mutual prime integers and &2a stands for the 
Brillouin phase. Rational values of the commensurability parameter, result 
in a finite number of energy bands. When it become an irrational number, 
the energy bands becomes infinite, the spectrum depending completely on the 
potential strength played by the anisothropy parameter. 

So, for A > 1 we have localized solutions, for A < 1 we have extended 
wave functions of a contiguous spectrum and for A = 1 we have a singular 
continuous spectrum. 

The Harper equation which originates from the minimal substitution on 
the energy dispersion law, can be rewritten in terms of transfer matrices [30] 
as 



where the transfer matrix [31] is given by 



> I V" "o 1 ) 



while the potential is periodic V n (x) = V n (x + 1) 



3 Numerical studies 

In order to investigate the energy bands when the system passes from a 
conductive state to an insulator one, we chose a particular class of non-analytic 
potentials, namely the Kohmoto potential [32J, proposed independently by 
Ostlund [33], because it preserves discontinuities and takes only two values 

f V A = 1 for - 7 < x < - 7 3 
V n = a = \ 5 
[ Vb = — 1 for — 7^ < x < 7^ 

where 7 = (\/5 — l)/2 represents the golden mean ratio. The Kohmoto model 
was introduced at the beginning as a simple toy model [Ml 131)] . with the 
purpose to explain the metal-insulator transition. By considering only two 
kind of possible values of V n , the Fibonacci sequence can readily established. 
Accordingly, the potential becomes [36] 

V n = V B + V A ([(n + 1)(3\ - [n(3\ ) (6) 
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where \n(3\ represents the greatest integer lower than nj3. Following the Nau- 
mis and Rodrigues work and using the identity n(3 = \ n{3\ +{n/3} the potential 
can be rewritten as 



Vn = VaP + Vb(1 - 0) + (V A - V B ){{n(5} - {(n + 1)(3}) 



(7) 



where {n/3} denotes the decimal part. 

In our numerical calculation the quantity (Va — Vb) in ([7]) known as the 
strength of the quasi-periodicity is taken as (Va — Vb) = — 1. The next step 
is to introduce this potential into the Harper equation, such as 



because the last one has many convenient symmetries, which can facilitate the 
numerical analysis. In this way the metal-insulator transition can be studied 
on the two-dimensional lattice, as a splitting effect of the energy bands, corre- 
sponding to the electrons and holes concerning energy bands. The numerical 
results are presented in Fig.l for several values of n as follows: 

a) The energy bands are continuous. The conductive medium exhibits 
mixed states of both electrons and holes if n = 1. 

b) One deals with splittings in the energy spectrum, supporting the cosine 
shape of the dispersion law of the energy if n = 1.1. 

c) The energy gap tends to increase releasing the Peierls gap if n = 1.2. 

d) More and more splittings occur in the inner bands distribution releas- 
ing incrementally the Peierls gap. The system is still a conductive medium if 



e) One deals with the Aubry duality point, which corresponds to the 
transition state from conduction to the insulator one. The Fermi level is 
located at the middle of the Peierls gap if n = 2. 

/) Beyond the critical point the energy gaps tend to increase. By now 
the system become more and more an insulator for n = 3. 

g) The lower energy bands tend to decrease in an accentuated chaotic 
manner, while the top bands tend to gather together around the zero value if 
n = 4. 

h) The separations between the energy bands become more pronounced, 
while the lower energy bands become a cloudy like configuration of points and 
short range lines if n = 7. 

i) Separations between energy bands are clear. Since an environment 
can't be a conductive medium with only one carriers, determine the insulator 
states. In this case state, the Peierls gap become to large between the energy 
of the negative carriers and the positive ones, to allow a combination. The 





n = 1.5. 
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conductive environment remain only with one kind of cherge carriers, which 
leads to the suppression of current flow if n = 100. 



Fig.l The energy spectrum versus (5 showing 
the band splittings for several n-values. 



By virtue of this transition the system exhibits less symmetric states, 
but with a lower-energy configuration. The energy bands splits, due to the 
periodic potential, which induce a small gap by releasing and expanding the 
Peierls distortion. The band splittings take place in such a way that the Fermi 
energy is always placed in the middle of the gap. The higher energy bands 
corresponding to the electrons get gathered araund the zero energy, while the 
negative ones corresponding to the holes lower their energy. This interpreta- 
tion of the energy band splittings is an atempt to investigate numerically the 
way in which systems undergo transitions from conductive to insulator states. 
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4 Conclusions 

In this paper we analyzed numerically the way in which the energy bands 
get splitting, when the system goes from conductive to the insulator state, 
i.e. to a less symmetric but lower-energy configuration, passing through the 
transition point. The energy bands become divided by increasing small gaps 
in such a way that the Fermi energy will be placed always in the middle of the 
gap. As result of n continuous variation, the system passes from a symmetric 
state to symmetry broken ones. The higher energy bands become locked, while 
the lower energy bands tend to become unlocked. 
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